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Abstract 

Boundary feedback stabilisation of linear port-Hamiltonian systems on an 
interval is considered. Generation and stability results already known for lin¬ 
ear feedback are extended to nonlinear dissipative feedback, both to static 
feedback control and dynamic control via an (exponentially stabilising) non¬ 
linear controller. A design method for nonlinear controllers of linear port- 
Hamiltonian systems is introduced. As a special case the Euler-Bernoulli 
beam is considered. 
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1 Introduction 

Within the last years there has been done a lot of research on the stability and 
stabilisation of wave and beam equations. Sufficient conditions which are easy to 
check for asymptotic or even exponential stabilisation of systems of the abstract 
port-Hamiltonian form 

^(t,C) = XPfc^^(t,C), t>0, ce(0,i) (1) 

where (K = R or C) via suitable dissipative linear boundary condi¬ 

tions have been given. Here, for the case = 1 (i.e. in particular wave equation, 
transport equation and Timoshenko beam equation) we should mention [19], the 
PhD thesis [18] and the monograph [9]. Also note the more recent article [7] for 
the case N = 1 and TL = I. More recently, investigations have been done in three 
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generalising directions. First, stabilisation using a dynamic controller ([14], [2]), sec¬ 
ondly considering the case TV > 2 ([2]) which includes the Schrodinger equation and 
the Euler-Bernoulli beam equation and last but not least nonlinear feedback ([11], 
[17]). This article should be seen as a continuation of [2] in which we extend the re¬ 
sults given therein for the linear feedback case to the situation of nonlinear feedback. 
Since in [2] the main tool for investigating stability were the Arendt-Batty-Lyubich- 
Vu Theorem and the Gearhart-Greiner-Priiss Theorem which both only hold for the 
case of linear evolution equations we have to develop alternative tools to attack the 
nonlinear feedback problem. Note that the infinite-dimensional system itself re¬ 
mains linear, so we do not touch the topic of nonlinear port-Hamiltonian systems 
(in the strict sense). Still, nonlinear feedback urges us to consider nonlinear contrac¬ 
tion semigroups (see, e.g. [13] and [16]) instead of linear semigroups, following ideas 
similar to those in [17] for the generation results and then exploiting ideas already 
used in [6] (in a linear scenario) for stability properties. We should also mention 
that the approach of [18] (there TV = 1 and linear feedback) may be used to obtain 
stability results for both the static and dynamic scenario (also see [11]). However, 
most likely this method is restricted to the case TV = 1. Also we stress that the 
interest in nonlinear feedback is motivated by applications where sometimes (due 
to physical or technical restrictions) it is not possible to implement perfectly linear 
controllers. Somehow the results of this article therefore show that to some extend 
nonlinear perturbations (from the perfectly linear case) do not harm the stabilisa¬ 
tion properties. Also note that in some cases the (usually finite dimensional) control 
systems considered here actually consist of both a finite dimensional controller and 
a finite dimensional control target which are connected mechanically via a beam 
modelled by a infinite-dimensional port-Hamiltonian system, e.g. wave equation, 
Timoshenko beam or Euler-Bernoulli beam, see e.g. [11]. 

This article is structured as follows. In Section 2 we recall some background on non¬ 
linear contraction semigroups and m-dissipative operators on Hilbert spaces, as well 
as the Komura-Kato Theorem which as generation theorem for nonlinear contrac¬ 
tion semigroups may be seen as the nonlinear analogon to the Lumer-Phillips The¬ 
orem. Then in Section 3 we recall and stress some properties of port-Hamiltonian 
systems in impedance passive boundary control and observation formulation. These 
observations together with results of Section 2 then give the generation result The¬ 
orem 4.3 (cf. Theorem 5.4 in [17]) which is restricted to static boundary conditions. 
In fact, we prove that the port-Hamiltonian operator A associated to the evolution 
equation (1) with suitable dissipative boundary conditions generates a (nonlinear) 
contraction semigroup (5(t))t>o on D{A){ = X ■.= A^(0,1;K'*) under suitable as¬ 
sumptions). From there we first consider the case TV = 1 and static feedback in 
Section 5 which serves as a introductory model case for the more general results 
later on. We are particularly interested in exponential stability of the nonlinear 
semigroup {S{t))t>o (in other words, global exponential stability of the equilibrium 
cc = 0), i.e. we ask whether there are constants TH > 1 and w < 0 such that for all 
X G D(A)(= X in this section) one has 

||S'(t)x||^ < TWe"^* ||x||^ , t>0. 

We show that under similar conditions as in the linear feedback case, exponential 
stability can be ensured. As example we consider boundary stabilisation of the 
wave equation. In Section 6 we leave the static feedback setting and consider a 
dissipative interconnection with a nonlinear dynamic controller. The main result 
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of that section is the generalisation of Theorem 4.3 to its dynamic counterpart 
Theorem 6.2. Then the subsequent sections are devoted to a) transferring the 
results of Section 5 to the dynamic scenario and b) generalising these results (both 
the static and the dynamic feedback cases) to port-Hamiltonian systems of order 
N = 2 where we also have a look on the Euler-Bernoulli Beam as a special case 
(where due to structural assumptions the dissipation assumptions for stabilisation 
are less restrictive). 


2 Some Background on Contraction Semigroups 

Before actually starting with the investigation of port-Hamiltonian systems we first 
recall some well-known concepts and results on semigroup theory. Since we only 
consider dissipative systems here, we restrict ourselves to the contractive case. For 
the general theory, see [13], and for the linear case see, e.g. the monograph [8]. 
Throughout we use the following definition. 

Definition 2.1 (Semigroup). Let X be a Banach space and Xq G X a closed subset. 
A family {S{t))t>o of mappings S{t) : Xq —>■ Xq {t > 0) is called semigroup if it 
satisfies the properties 

1. S'(O) = Ixo 7 ths identity map on Xq, and 

2. S{t -I- s) = S{t)S{s) for all s,t > 0. 

We speak of a strongly continuous (abbr.: s.c.) (nonlinear) semigroup (or, dynam¬ 
ical system) if S(t) € C{Xq;Xq) (t > 0) and for all x € Xq the map t i—>■ S(t)x 
is continuous on R_|_ := [0, oo). If additionally Xq = X and all maps S{t) (t > 0) 
are linear, i.e. S{t) € B{X), then we speak of a strongly continuous semigroup (of 
linear operators), or Co-semigroup. A (linear or nonlinear) semigroup {S{t))t>Q is 
called contractive, if all maps S(t) (t > 0) are contractions, i.e. 

||S'(t)a; — 5'(t)a:'||^ < ||a: — , x,x' € X, t>Q. 

Below we restrict ourselves to the Hilbert space case. We are going to state the 
nonlinear version of the Lumer-Phillips Theorem, i.e. the Komura-Kato Theorem, 
and therefore recall the concepts of dissipative (resp. monotone) operators. For 
details see, e.g. Chapter IV in [16]. 

Definition 2.2. Let A : X ^ U{X) := {B C X} be a map. We write D{A) for its 
domain 

D{A) := {x € X : A{x) ^ 0} 

and also write A : D{A) C X =4 X. If A{x) = {j/x} for all x £ D{A) we call A an 
operator and write 


Ax := px, for the unique yx G A(x). 

Otherwise we say that A is multi-valued. 

We use the notation H -|- i? for the sum of two maps A : D{A) C X ^ X and 
B : D{B) C X =4 X as follows 


{A B){x) = {yi y 2 £ X : yi £ A{x), y 2 £ B{x)}. 


PH SYSTEMS WITH NONLINEAR BOUNDARY FEEDBACK 


4 


Note that D{A + B) = D{A) n D{B) and for the particular case where B is an 
operator {A + B){x) = {yi + Bx : yi G A{x)}. 

Definition 2.3. Let X be a Hilbert space and A : D{A) Q X ^ X. We say that 
A is dissipative (and —A monotone (or, accretive)), if for all x,x' G D{A) and 
y G A{x), y' G A{x') one has 

Re {y-y' \ x- x')^ < 0. 


If additionally 


{y G X :3x G D{A) y G {I — ^)(x)} =: ran {I — A) = X 


then A (resp. —A) is called m-dissipative (resp. m-monotone (or, m-accretive)). 
If A : D{A) G X ^ X is dissipative and has no proper dissipative extension, i.e. if 
B : D{B) G X X with D{A) C D[B) and A{x) C Bix) for all x G D{A), then 
A = B, then we call A maximal dissipative. 

Remark 2.4. Let A : D{A) X be an m-dissipative map on a Hilbert space X, 
then for all x G D{A) the set A{x) is convex and thus there is a unique z G A(x) 
with minimal norm. This defines the minimal section A^ of A: 

A°x:=z, ||z||= inf bll, D{A°) = D{A). 

y£A{x) 

Moreover for all x G X and X G the element y G D{A) such that 

xG {XI - A){y) 

is uniquely determined, thus we may write y = {XI — A)~^x. In particular every 
m-dissipative operator is maximal dissipative. 

Lemma 2.5. If A : D{A) X is m-dissipative and B : X ^ X is dissipative and 
Lipschitz continuous, then also A-\- B : D{A) X is m-dissipative. 


Proof. See Lemma IV.2.1 in [16]. □ 

For clarity, let us also mention Minty’s Theorem. 

Theorem 2.6 (Minty). On a Hilbert space X the m-dissipative operators are exactly 
the maximal dissipative operators. 


Proof. Combine Lemma 2.2.12(iii) and Corollary 3.2.27 in [13]. □ 

As in the linear (Cp-semigroup) case, m-dissipative operators are closely related to 
the generators of contraction semigroups. 

Definition 2.7. Let {S{t))t>o be a nonlinear strongly continuous contraction semi¬ 
group on X. Set 

D := {x G X : S{-)x G Lip(]R+; V)} . 

We define the (infinitesimal) generator of the s.c. contraction semigroup {S{t))t>o 
as 

. / N .. S{t)x — x , r S{t)x — X , 

Aq{x) := lim-, D{Aq) := Ix G X : lim-G X exists} 

t\,o t t\o t 

and the (g)-operator A : D{A) C V X as the maximal dissipative extension of 
Aq with D{A) G D. 
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Remark 2.8. By Zorn’s Lemma every dissipative operator has a maximal dissipa¬ 
tive extension (see Lemma 2.2.12(ii) in [13]). Hence the (g)-operator always exists. 
Also note that the infinitesimal operator Aq (or the (g)-operator A) uniquely deter¬ 
mines the s.c. contraction semigroup (see Corollary 3.fill in [13]). 

Theorem 2.9 (Komura-Kato). Let A : D{A) C X X be a (possibly multi¬ 
valued) map on a Hilbert space X. If A is m-dissipative, then it generates a 

_ ^ 

nonlinear strongly continuous contraction semigroup {S(t))t>o on X_ := D{A) . 
More precisely, for each xq € D(A) there is a unique absolutely continuous solution 
X) of the abstract nonlinear Cauchy problem 

G A{x(t)), t > 0 

x(0)=xo. (2) 


Also ||^a;||^ (R+-x) — tho function 

t > 0 and the right-derivative ^ one has 

^x(i) = A°x{t), t>0. 


is decreasing and for every 


Proof. See Proposition IV.3.1 in [16]. □ 

Remark 2.10. If A is m-dissipative and 0 G ^(0), then S(t){Q) = 0 for all t > 0. 
Consequently in this case 

||5'(t)x||;^ < Ijxll;^, t>0. 


3 Impedance Passive Port-Hamiltonian Systems 

In this section we lay the foundations for the generation theorems later on as we 
introduce port-Hamiltonian systems and boundary control and observation systems. 
We assume impedance passivity (as boundary control and observation system) and 
observe that the transfer function exists on K)}" and its symmetric parts are coercive 
as linear operators on 

Definition 3.1 (Port-Hamiltonian System). Let TV G N and Pk G {k = 

0,1,..., N) with Pf = Pk (fc = 1,..., N) and the symmetric part Re Pq := 

< 0 of the matrix Pq be negative semi-definite as well as H. £ Too(0,1)'^^'^ 
and consider the Hilbert space X = ^2(0,1)"^. Further let Wb,Wc G he 

matrices such that G invertible. 

1. If Pn is invertible and % is coercive as multiplication operator on X, i.e. 
there is mg > 0 such that 

z*H{C)z > mo \z\^ , z G a.e. ( G (0,1) 

then the operator 

mx = J2Pk-^i^x) 

k^O ^ 

D(2l) = {a; G L2(0, l)"^ : Ux G R^(0,1)^} 
is called (maximal) port-Hamiltonian operator. 
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2. For a port-Hamiltonian operator 21 the boundary port-variables fd,'Hx G 
(^boundary flow^ and eg^-ux S ('boundary effort^ are defined as 


/ {nx)(i) \ 


[ fd,'Hx 

V ^d.Hx 


= Re. 


{Hxm 




where Rext = ^ ^ 2 Ndx 2 Nd (invertible) matrix 


Q 


Pi P2 . Pn ■ 

—P2 —P3 ■■■ —Pn 0 


0 ••• 0 0. 


3 . For a port-Hamiltonian operator 21 we define the input map 05 : i0(2l) C X —>■ 
and the output map £ : 10 (2t) C X —>■ via 

and call & = (21, 05, £) a port-Hamiltonian system to which we associate the 
following abstract boundary control and observation problem 

—a;(i) = 2la;(t), a;(0) = xq 

u{t) = iBx{t) 

y{t) = £x(i), t>0. 

From here on we assume that © = (21, 05, £) is a port-Hamiltonian system in the 
sense of Definition 3.1. Further we always assume that the system © = (21, 05, £) is 
impedance passive. 

Assumption 3.2. The port-Hamiltonian system © = (21, 05, £) is impedance pas¬ 
sive, i.e. 

Re (2ta; | x)-^ < Re {3Bx \ €x)^Nd , x G D(2l) 
where we take X = L2{0,1)‘^ to be equipped with the inner product (• | ■)^ := 

Remark 3.3. Assumption 3.2 already implies that the operators 2l|ker® ond 2t|kerC 
are dissipative on X. Hence by Theorem 7.2.) in [9] these operators generate s.c. 
contraction semigroups (of linear operators), respectively. In fact, below we prove a 
nonlinear version of this result, see Theorem ).3. 

One may even say more, namely © = (21, IB, £) also is a boundary control and 
observation System. 

Proposition 3.4. The impedance passive port-Hamiltonian system © = (21, 05, £) 
is a boundary control and observation system, i.e. 
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1. 21 : £>(21) C X ^ X, «8 : £>(21) C X ^ U, £ : £>(2l) C X ^ Y are linear 
operators where X, U and Y are Hilbert spaces. (Here: U = Y = 

2. 2l|ker® generates a Co-semigroup on X. (Here: The semigroup is even con¬ 
tractive.) 

3. There is a right-inverse B € B(U,X) of^ such that 

ran B C £)(2t), 2t£ G B{U,X), »£ = I. 


4 . € is bounded from ker® to Y where ker® is equipped with the graph norm 
o/2l|ker<8- (Here: £ is also bounded from D(2t) to Y where £>(2t) is equipped 
with the graph norm of the (here: closed) operator ill.) 


Proof. See Theorem 4.4 in [12]. 


□ 


Remark 3.5. If Wb 


Wb,i 

^B,2 


for Wb,i G g ^{Nd-k)x2Nd 


k G {0,1,..., Nd} (and an analogue decomposition of Wc, ® and €, respectively) 
and one considers 


2li := 2lb(ai), D{^i) ■= {x G £(2t) ; 3&2X = 0}, 

then (quite naturally) also (2ti, ®i, £ 1 ) is a boundary control and observation system 
(where Ui = Yi = since this only means that we fix some components of the 
input u to be zero. Hence without loss of generality we always assume k = Nd. 
(Here we use the notations K° := {0} and := S(K”,K™) also in the cases 

where n or m equals zero.) 

Let us also introduce the concept of a transfer function which is closely related 
to the Laplace transform of the semigroup (T{t))t>o generated by A. (For more 
background on transfer functions for port-Hamiltonian systems we refer to Chapter 
12 in [9].) 

Definition 3.6 (Transfer function). Consider the abstract boundary control and 
observation problem 


^x{t) = 2la:(t) 

u(t) = ®x(f) 

y{t) = £x{t), t>0 

where © = (2t, ®, £) is a boundary control and observation system and let A G K. 
We write A G D{G) if there is G(A) G B{U,Y) such that for all u G U there is a 
unique solution of 


Xx = 21a; 
u = ®a; 

y=£x (3) 


where x G D(2l) and y G Y is given by y = G(X)u. 
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For impedance passive port-Hamiltonian systems one has Kj]" :={AsK:ReA> 
0} C D{G). 

Lemma 3.7. Let & = (2t, *8,£) be an impedance passive port-Hamiltonian system 
with H = I. Then Kj]" C D{G) and Re G'(A) > 0 for all A S Kj, i.e. for all A € 
there is m\ > 0 such that 

Re (z I G{X)z)^ > mx jzj^ , z€U = 

More precisely, for every A € Kj there are operators $(A) G B(N), 'lt(A) G B(U,X) 
and F{X) G B{X, Y) such that for all f G X and u G U there is a unique solution 
of the problem 


(A -%)x = f 
u = iBx 
y = <Lx 


which is given by 


X = $(A)/ + 'I'(A)'U 
y = F{X)f + G{X)u. 


Remark 3.8. The restriction % = I is not necessary. In fact, for any impedance 
passive boundary control and observation system © = (21, 25, £) (on Hilbert spaces X 
and U = Y) and P G B{X) any coercive operator on X, also 6p = (2tP, *8P, £P) 
is an impedance passive boundary control and observation system (on Xp = X 
equipped with (• | ■)xp ■= (’ I P’)x) transfer function exists on ]K(}" for 

(21, 25,£) if and only if it exists on for 6p = (2tP, *8P, £P) (the situation is 
similar for $, and F as in Lemma 3.7). 


Proof of Lemma 3.7. Let A G K^J", u G and f G X be given. First, observe 
that the equation 

(A - 2t)x = / 

has the general solution x = for ^ := (x,x',... and ^(C) = e‘“^^^(0) + 

9/(0 where 

r 0 1 0 ■■■ on 


Pa 


0 


0 1 


6 . 

-Pn^Pi - 


0 

0 1 
- -p;^Wn-i_ 


and qfiC) = 
expressed as 



ds. Writing Ex = e^^ input and output may be 


U = WsRecct [5^00) + WpRext , 

y = WcRe.t [ t ] OO) + WcRe.t . 


Since the system (2t, 25,£) is impedance passive both the matrices WpRext 
and WcRext are invertible since otherwise (choosing / = 0 and ^(0) in the 
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kernel of one of these matrices) A G Kq fl cr(2t|ker<8) or A G ]&)}■ fl cr(2l|kerir), hr 
contradiction to 2l|ker<8 and 2t|kerC being dissipative. As a result, for any given 
u G and f G X there is a unique solution {x,y) G £*(2t) x and clearly 


the map (/, tt) {x,y) =: [1] is linear and bounded. By the same 

reasoning one hnds (for / = 0 hxed) the inverse map G(A)“^ : y u, so that G(A) 
is bijective. Further we have for all u G \ {0} and the corresponding solution 
{x,y) G d?(2l) X of (3) that 

Re {u I G{X)u)^Nd > Re (Qlx | x)p^ = Re (Ax | x)p^ = Re A ||x||^^ > 0 
so that in fact the symmetric part Re G(A) > 0 is strictly positive definite. □ 


4 A Generation Theorem for Static Feedback 

To begin with we prove the generation theorem for port-Hamiltonian systems with 
nonlinear dissipative boundary conditions. We use a strategy very similar to the 
linear case, in fact the main differences are twofold: On the one hand we use the 
nonlinear generalisation of the Lumer-Phillips Theorem, namely the Komura-Kato 
Theorem, so that again the proof of the generation result reduces to showing (besides 
dissipativity which is an assumption) a range condition. In the linear case it proved 
convenient (however, not necessary) to reduce the generation theorem to the special 
case where TL = I. As an additional hurdle, the relevant Lemma 7.2.3 in [9] has to 
be formulated in a nonlinear version, and regarding the proof one should replace 
any reasoning with the (linear) adjoint by an argument which is suitable for the 
nonlinear situation. In fact, this is 

Lemma 4.1. Let X he a Hilbert space and A : D{A) Q X X be a dissipative, 
possibly nonlinear and multivalued, map. Further assume that P G B{X) is strictly 
coercive. Then if A — I is surjective, so is AP — I. 

Remark 4.2. Note that this a very special and simple case of Theorem 2 in [3]. 
Since the proof of Lemma j.l is quite elementary we give it nevertheless. 

Proof of Lemma 4.1. As a first step, assume that ||P —/|| < then from 
Neumann’s series we conclude that ||P~^|| < i_||p_j|| < = 2 so that 

||P-/||||P-1|| =:pG (0,1). 

We shall show that for any given f G X there is x G D{AP) such that 

{AP-I){x)^f 

which is equivalent to solving the problem 

[AP - P){x) 3 f + {I - P)x, 


or, since {A — I) ^ exists, 

X = 4)/(x) := P~^{A — I)~^ (/ + (/ — P)x). 
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We show that : X —X is a strict contraction and therefore admits a unique 
fixed point Xf =: {AP — I)~^f. In fact, we have 

< ||P-1|| ||(A - /)-!(/ + (J - P)x) -{A- I)-\f + (/ - P)x')\\ 

< \\P-^\\\\if+{I-P)x)-if + iI-P)x')\\ 

< ll-P"! =p\\x-x'\\ 

where we used Corollary 1.3(b) in [16] in the second step. Therefore $/ is a strict 
contraction and the Contraction Principle gives a unique solution Xf =: (AP — 
I)~^f. In the second step we remove the restriction on P. Namely it is easy to see 
that there are a number n € N and a coercive operator Q = g }3{X) such that 
||/ — Q\\ <\ and P = Q". Note that ||/ — Q\\ = ||/ — Qjj^, for all the norms 



fc = 0,1,..., n — 1. 


Writing 


AP-I = - / 


the general case follows by induction using the spaces Xk := {X, IHIj.). □ 

It is an easy consequence of the preceding lemma that for any m-dissipative operator 
A : D{A) X ^ X and a strictly coercive operator P G BiX) also the operator 
AP with domain D{AP) = {x € X : Px G D{A)} is m-dissipative on X equipped 
with the equivalent inner product (• | •)p = (• | P-) (cf. Lemma 5.1 in [17]). Of 
course, in our particular situation P = H. is the Hamiltonian density multiplication 
operator. 


Theorem 4.3. Let & = (21, 25, £) be an impedance passive port-Hamiltonian sys¬ 
tem. Assume that cj) : is a (possibly multi-valued, nonlinear) m- 

monotone map. Then the (single-valued) operator 


^ = ^\d{A) 

D{A) = {a; G P(2l) : 25a; G -(/-(Ca;)} 


generates a s.c. contraction semigroup on X = ^2(0,1)“^ with the inner product 

Remark 4.4. Note that for the case N = 1 a characterisation of m-dissipative 
boundary conditions yielding an m-dissipative operator A has been given in Theorem 
5.4 of [17]. Also note the more general result Theorem 3.1 therein. 


Proof of Theorem 4.3. From Lemma 4.1 we know that it suffices to consider the 
case where H = I equals the identity. Also note that there is xq G D{A) ^ 0 which 
implies that xo-l-C'))°(0, l)'^ C D{A) is a dense subset of X. Clearly A is dissipative 
since for x,x G D{A) we have 

Re {A{x) — A{x) 1 x — x)^^ = Re (2t(a; — x) [a; — x)^^ 

< Re (25a; — 25a; ] €x — <I.x)^Nd < 0 
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using that 58a; G —iBx G —cj}{<Lx) and is monotone. It remains to show 
that ran {I — A) = X, i.e. for every f G X we have to find x G D{W) such that 

(/ - 2t)x = / 

58a; G —(j){<Lx). 

From Lemma 3.7 we know that all solutions of the first of these equations have the 
form a; = where 


e(C) =e«^^e(0) + 





ds 


and the problem thus reduces to hnding ^(0) G 'K^‘^ such that 


WBRe.t [f]^(0) G -(l>{WcRe.t[j]m) (4) 

where E = Ei = or, by Lemma 3.7, finding y G such that 
u = G(l)-'y-G(l)-iF(l)/G-<^(y), 

i.e. (G(l)“^ + (/>)(y) 9 G(1)“^F'(1)/. Since </) is m-monotone and Re G(l)“^ is 
coercive by Lemma 3.7, also </> + G(l)“^ — el is m-monotone by Lemma 2.5 for 
some small e > 0. We conclude that there is a (unique) y G such that for 
u := G{l)~^y one has u G —(j^iy) and since the matrix WsRext [f ] is invertible 
it follows that there is a (unique) ^(0) G such that (4) holds, i.e. there is 
a (unique) x G D{A) with f G {I — ^)(a;). Now the assertion follows from the 
Komura-Kato Theorem. □ 


5 Stabilisation of First Order Systems 

For the moment we additionally assume that N = 1, i.e. ^ = Pi{'H-y + Po{'H-) on 
£1(21) = {xG £ 2 ( 0 , l)'^ : Hx G H^iO, 1)'^}. Also we assume that H G W^lo, 
is Lipschitz continuous. 

We aim to prove the following uniform exponential stability result. 

Theorem 5.1. Let © = (21, 58, £) be an impedance passive port-Hamiltonian system 
and (j) : ^ an m-monotone map with 0 G (/>(0). For the nonlinear operator 

A := 21|d(a), D(A) := {x G £>(2l) : *8x G -(j^i^x)} 

assume that there is k > 0 such that 

Re {Ax I x)-f^ < —k{x*TLx){1), x G D(A). 

Then A generates a s.c. contraction semigroup {S{t))t>o with globally exponentially 
stable equilibrium 0, i.e. there are M >1 and w < 0 such that 


||S'(t)a;||„ < ||x||„ , 


X G X, t>0. 
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Remark 5.2. If (j) e B{U) = is linear this is Theorem HL2 in [19] which uses 
a “sideways energy estimate” (Lemma HI.l in [19]) in the spirit of [5[. Actually a 
first result like this may already be found as Theorem 3 in [15] where H is assumed 
to be smooth. An alternative proof of Theorem HI.2 in [19] via Gearhart’s Theorem 
can be found as Proposition 2.12 in [2], but clearly the latter technique is restricted 
to the linear situation. In fact, we will use a technique which for the linear case 
yields a third proof of the theorem. Namely we use an idea which was used in [6] to 
prove exponential stability for a chain of linear Euler-Bernoulli beams with (linear) 
dissipative linkage. 

As preparation we state the following auxiliary result. 

Lemma 5.3. Let a > 0 and / 3,7 > 0 be given. Then there is rj G C°“([0,1];M) with 
r]{0) = 0 and rj' > 0 such that 

arj'iO - PviO >1, Ce[0,1]. (5) 

Proof. Scaling rj by the factor ^ it is enough to consider the case 7 = 1 . We make 
the ansatz 77 (C) = 6 ^^ ~ 1 foi' A > 0 which we are going to specify. Then equation 
(5) is equivalent to 

(oA — /3)e'^^ > 7 — /3 (C G [0,1]) <t7 oA > 7 

A>^. 

a 


□ 

Also the following fact (which can be derived via integration by parts) will prove 
quite useful in the computations below. 

Lemma 5.4. Let Q G W^(0,1)'^^‘^ be a function of self-adjoint operators and 
uGH^{ 0,1)^. Then 

Re [u' I Qu)^^ = (u I Q'u)^^ + i [m(C)*Q(C)m(C)]o • 

Proof of Theorem 5.1. Existence of the s.c. contraction semigroup {S{t))t>o 
follows from Theorem 4.3. Stability: On X = ^ 2 ( 0 , 1 )"^ define the quadratic func¬ 
tional 

q{x) := {x\r]PC^x)^^, xGX 

where 77 G C'^([0,1]; R) is a function with 77 ( 0 ) = 0 and 77' > 0 which we choose at a 
later point. Let xq G D(A) be arbitrary and denote by x(t,-) = x(t) = S(t)xo the 
solution for the initial value xq, so that x G Wt^(R+; A), hence q{x) G MOo(R-i-;R) 
and we calculate (using Lemma 5.4) for a.e. t > 0 

j^q{x(t)) = 2Re {PC^xif) I vx(t))^^ 

= 2Re {{(Hx{t))' P([^PoiTdx{t)) \ pxft)) 

= -{nx{t) I YriifH-^)' -2riRe{n-~^PC^PQ))'Hx{t))^^ 

+ [nx{tx)*{pn-^){C)nx{t,o\l. 
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Thus for $(f) := t ||a::(f)||^ + q{x(t)) we have 

= I {ii-ri^)n-^ + rj[-{n-^y + 2Re {n-^pr^Po)])nx{t))^^ 

+ ri{l)x{t, l)*'H{l)x{t, 1) + 2fRe {Ax{t) \ x(t))^ 

< {Hx{t) I ((l-77')'H-^+r?[-CH-^)' + 2Re {R-^Py^P q)]) Hx{t)) 

+ (??(1) — 2tK)x*{t, l)H{l)x{t, 1). 


So far, we did not specify our choice of 77 , so we may choose 77 in such a way that 


(l-77'(C))?^-i(C)+r7(C) [(H-i)'(C)+2Re {n~\OPC^Po)] 

< (Mo - v'iOmo + viO [Ml + 2 M 3 ]) / < 0, a.e. C £ (0,1) 


where 


mol < n~\c) < Mol, (n-^YiC) < mr 
-M ol < Re {n-\OPy^Po) < Mol 

for a.e. C £ (0,1) and we applied Lemma 5.3. Then for t> to '■= we have 

lm<o 

and thus $ decreases on (to, 00 ). (Note that the choice of 77 and to is independent 
of the initial value xq £ D(A).) Moreover, since \q{y)\ < c|| 7 /||.j^ for some c > 0 and 
all y £ X we obtain for t > to the estimate 

t\\x(t)\\i<m+c\\xit)f 


H 


< ^(to) + c\\x(t)\\l^ 


and hence for t > max{to, c} 


Mt)\\n< 


2 / ^(io) ^ to + c 


< 

t — c t — c 


<to)\\ 


H 


< ^ Ikolll 


and from the density of D(A) in X we conclude for t > max{to, c} 


fo + c 


t^+oo 


t — C 


\x 


> 0, X G X. 


As in the linear case, this property easily implies uniform exponential energy decay. 

□ 


Remark 5.5. For q as in the proof of Theorem 5.1 the following holds. For every 
solution X G W(^(IR.+ ;X) fl Loo(K+; d?(2t)) of x = 2lx one has q(x) G W,^(M+;M) 
with 

\\^(y)\\n + ^ 9 (a;( 0 ) <c\(nx)(t,l)\^ , a.e. t > 0 . 

We come back to this property in the eontext of dynamic controllers. 
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Remark 5.6. An alternative proof can be established via the “Sideways Energy 
Estimate” of Lemma 9.1.2 in [9]. In fact, the proof of Theorem 9.1.3 in [9] almost 
literally carries over to the nonlinear situation. 


Let us also state an asymptotic stability result which follows from Theorem 5.1, the 
contraction property and the following interpolation inequality. 

Lemma 5.7. Let 1 < j < fV £ N fee natural numbers, then there is a constant 
C > 0 such that for every f £ H^ {0, 1) one has 


fU) 


C[0.1] 


< C 


(0.1) 


+ 


fiN) 


i2(0.1) 


27 + 1 
2iV 


I 2N 
1 ^ 2 ( 0 , 1 ) 


Proof. Combine Lemmas 4.10 and 4.12 in [1] to get a one-dimensional version of 
Theorem 4.14 therein. □ 

Remark 5.8. Actually Lemma 5.5 together with Theorem IV.1.1 in [10] implies 
that 21 is a closed operator. 

Corollary 5.9. Under the assumptions of Theorem 5.1, but with the less restrictive 
condition 


Re {Ax I a;).^ < —k{x*'Hx){1), x £ D{A) : |*8a;| , \^x\ < p 

for some p > 0, the semigroup {S{t))t>o is asymptotically stable, i.e. for all x G X 
one has S{t)x 0. 

Proof. First step: Take any a;Q £ D{A) and set x{t) := S{t)xo. Then ||a;(t)||.j^ 
and ||2lx(fe)||.^ are bounded by ||a;o||.^ and ||2ta;o||^ for a.e. t > 0, respectively. 

Hence also |*Ba;(t)|, |£a;(t)| < c||a;o|| 2 i := cyj ||a;o||^ -b ||2ta;o||^ for a.e. fe > 0, so 
that x{t) = S{t)xo where S{t) is the s.c. contraction semigroup corresponding to 
A := where D{A) = {x G D{W) : *8a: £ —(j){€x)} for some m-dissipative 

]) : D{4>) C such that (^ = (/> for \z\ < c ll^olla and such that 

Re (^Ax I x^ < —k{x*'Hx){l), x G D(A). 

Consequently x(t) = S{t)xo 0 due to Theorem 5.1. 

Second step: Let x G X he arbitrary. We have to prove that for all e > 0 there 
is Tg > 0 such that ||S'(fe)a;|| < e for all t > T^. For this end take any x^ £ D{A) 
such that III — Xe II < §. Then by the first step for > 0 sufficiently large one has 
||5'(t)a;E|| < I (t > Te), so that 

||S'(t)x|| < ||S'(<)a; - 5'(fe)a;E|| -f ||5(t)a::E|| 

< Ik-+ ll-S'(0aie|| < e, t>T^. 


□ 


Example 5.10 (Wave Equation). Consider the one-dimensional wave equation 

puJtt{t, C) - {EIu;^)cit, 0=0, c e (0,1), t > 0 
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where EI,p G Loo(0,1) are uniformly positive, i.e. also El ^ G ioo(0,l). At 
the left end we assume conservative or dissipative boundary conditions of the form 


wt(t,0) = 0 or {Elu}(;){t,0) e f{ujt{t,0)), t>0 


where / : K. =| M is maximal monotone and /(O) 9 0, e.g. f could be single¬ 
valued, continuous and non decreasing with /(O) = 0, in particular the case / = 0 
(Neumann-boundary condition) is allowed. We further assume that on the right end 
a (monotone) damper is attached to the system, so that the boundary condition is 
given by 

{EIu)cO{tA) e -g{uJt{t,l)) 

where again g : R K is maximal monotone with g(0) ^ 0. Of course, the choice 
f = g = 0 would lead to Neumann-boundary conditions on both sides for which the 
system is known to be energy-preserving, in particular not strongly stable. Here as 
usual the energy is given by 

m := Pio 01 "+Eiio \u^^{t, 01 " dc 

Jo 

In fact the example fits into our port-Hamiltonian setting when we choose x = 
{puJt, ujz), H — disLg{p~^, El) and Pi = (^), Pq = 0 and K = R. If we define 




€x 


-{ElivOm -{nx)2{0) \ 

(P/u;c)(l) )~[ inx)2{l) ) 

^t(o) \_f {nx),{0) \ 

Ml) ) - 1 , {Hx)^{l) ) 


then the system & = (2t, S, £) is impedance passive, since for the maximal operator 
21 one has 


(21a; I a;)„ = {nx)i{l)*{nx) 2 {l) - {Hx)x{Q)*{Hx) 2 {Q). 

(Note that in the Dirichlet case Wt(0) = 0 one has to exchange the first components 
of 18 and £ and then choose f = 0.^ The corresponding port-Hamiltonian operator 
A = 11 I\d{A) (with nonlinear boundary conditions) is dissipative then 

P(A) = {a;GL2(0,l;KO; Pa: G P^(0,1;«"), (Px) 2 (l) G-5((Pa;)i(l)), 

J(Pa;)i(0) = 0, (Dirichlet b.c.) or 

|^(Px)2(0) G /((Pa:)i(0)) (Neumann b.c.) 

and we have at least 

(Ax I a:)^ < -(Pa:)i(l) 5 °((Pa:)i(l)), x G D{A). 

Theorem 4-3 assures that A generates a nonlinear contraction semigroup on X = 
p2(0,l;l^") with inner product (■ | •)^. To have stability results we need stronger 
assumptions on the damper, i.e. on the map g. First assume (additionally to g being 
m-monotone) that there is n > 0 such that k~^ |x| < \z\ < k |x| for all x and 
z G g(x) (i.e. in particular g(0) = {0}j. Then we obtain the dissipativity condition 

(Ax I a;)„ < -;t|(Pa;)(l)|" , 


X G D(A) 
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where k := ^ miii{K, and so Theorem 5.1 ensures uniform exponential stability 
of the corresponding nonlinear semigroup. Secondly we assume that the condition 
\x\ < \z\ < k\x\ only holds for a; G R and z € g{x) whenever \x\ < p for 
some fixed p > 0. Then we only obtain asymptotic stability of all solutions x(t) = 
S(t)xo {t > 0). We refer to Example 3.3 in [4] for suffcient conditions leading to 
strong stability of the n-dimensional wave equation on a smooth, bounded domain 
D C R". 

Example 5.11 (Dissipation near 0). Let (() : R'^ =4 R'^ be a map such that 

1 . (j) is m-monotone 

2 . |j/| > c\z\ for some c > 0 and all z in some open ball Bo{p) around 0 and 

y e fiz). 

Also let Ip : R'^ ^ R'^ be any m-monotone map and consider 

Ax = Qlx, D(Ql) = {xe D(Ql) : 'Hx[0) G V'(A'Hx(0)), 'Hx(l) G -(/)(PiHx(l))} 

where % is a first order port-Hamiltonian operator. Then A generates an asymptot¬ 
ically stable nonlinear contraction semigroup {S{t))t>o on X = L2(0,l)‘^ equipped 
with the inner product (• | 

6 A Generation Result for Dynamic Controllers 

Let us again consider a port-Hamiltonian system 6 = (21,58, £) of arbitrary order 
iV G N which is impedance passive. However, instead of a static feedback *Ba; G 
—4){Cx) we now consider the feedback interconnection 58a; = —yc and Uc = €x with 
a nonlinear control system Ec. 

To motivate the subsequent definitions and assumptions let us first consider the 
case of a finite dimensional linear system Ec = {Ac, Be, Cc, Dc) given by 

d 

-^^Xc{fi) — AcXc{h) -\- BcUc{fi) 

yc{t) = CcXc{t) + DcUc{t), t > 0. (LC) 

This system is impedance passive if and only if the matrix 



is dissipative (and then m-dissipative since Xc x is finite dimensional and the 
map is linear). So much for the linear and finite dimensional case. More general 
we assume that we have a Hilbert space Xc as the controller state space (equipped 
with some inner product (• | ) and Me : Xc x =1 Xc x is a (possibly 

multi-valued) m-dissipative map for which its minimal section generates a nonlinear 
s.c. contraction semigroup on Xc x 
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Example 6.1. If we have that Me = () where Ac : D{Ac) C Xc Xc and 
—Dc : D{—Dc) = ^ are m-dissipative and the operators Be : > Xc 

and Ce-.Xe^ are assumed to be linear and adjoint to each other, i.e. input 

and output are collocated, then Me : D(AIc) = D(Ac) x C Xc x =4 

Xc X is m-dissipative. (Note that choosing Be = 0 leads to static feedback 

as investigated before where the nonlinear system is decoupled from the infinite¬ 
dimensional part.) 

We denote by IIx^ : Xc x —>■ Xc and : Xc x —>■ the canonical 

projections on Xc and respectively. Then we may define A : D{A) C XxXc A 
X X Xc as 

H-XcMcixc^x) 

D{A) = {(a;,Xc) G D{^) x Ilx„D{Mc) ■ fBa; G Il^NdMc{xc,€x)} 

and we have the following generation theorem for the interconnected system. We 
use the notation 2Lc '■= Ux^D{Mc). 

Theorem 6.2. The (possibly multi-valued) map A : D{A) C X x Xc =4 X x Xc 
is m-dissipative on the product space X x 2Lc> thus its minimal section generates a 
nonlinear s.c. contraction semigroup {S(t))t>o on X x X^. 

Proof. From the Komura-Kato Theorem we know that it is enough to show m- 
dissipativity. Again we may and will assume that H = I. First let us show that 
D{A) = X X X^. Take any {x,Xc) G X x X^. As a first step, let us additionally 
assume that Xc G IlxcD{Mc). Then there are Uc and pc G 'K^‘^ such that {xc, Uc) G 
D{AIc) and pc G H^NdMc{Xc,Uc). We need to find a sequence {xn)n>i C I?(2l) 
converging to x (in X) and such that = Uc and = Pc- For this take any 
xo G 11(21) such that Siq = Pc and Cxq = Uc. Since C)!°{0A)‘^ is dense in X 
there is a sequence (z„)„>i C C^{0AY converging to x — xq (in X). Note that 
then Xn := Xo + Zn X with Sxn = *Bxo G nK„^Mc(xc, Cxq) = nK„^Mc(xc, £x„) 
does the job. As a second step we allow that (x,Xc) merely lies in X x Xc. Then 
there is a sequence {xc,n)n>i C IlxcD{Mc) such that ||xc,n — Xc\\x ^ n- Further 
we know from the first step that there are sequences {xn,k,Xc,n,k)k>i C D[A) 
such that \\{Xn,k,Xc,n,k) — {Xn, Xc,n)\\xxx — I hence the diagonal sequence 
(xn,n,Xc,n,n)n>i C D(A) converges to (x,Xc). Therefore D(A) is dense in X x 
Xc- Secondly, A is dissipative since for all (x, Xc), (x, Xc) G D(A), (Six, Zc) G 
A(x,Xc), (2lx, Zc) G A(x,Xc) we have 

Re ((2tx,Zc) - (2lx, Zc) | (x,Xc) - (x,Xc))xxx„ 

= Re (2l(x — x) I X — x)y^ + Re (zc — Zc | Xc — Xc)x 
< Re (^(x - x) I £(x - Xc))-KNd + Re (zc - Zc | Xc - Xc)^^ 

= Re ((gjx) — ((|j) I id) ~ iez))Xa ® 

since ( ^ ) and G Me ( ) for the m-dissipative map Me. Finally, 

we show the range condition ran(/ — .4) = X x Xc. Let (f,fc) G X x Xc be 
arbitrary. We look for (x,Xc) G D{A) such that 



(x,Xc) - (/,/c) G A(x,Xc) 
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which may be equivalently expressed as {x,Xc) G -0(21) x IIx^D(Mc) being such 
that 

(I -Ql)x = f 

mx € njjiv.M, 

where from the first equality and Lemma 3.7 we get x = $(!)/ +and 
£x = E(l)f + G(l)^x. Since 0(1) is invertible it only remains to solve the problem 

(c(1?V.)-(g(1)4(1)/)'^"‘(S)- (•=) 

Since for some e > 0 small enough el — Re G(l)~^ is still dissipative, clearly 

A - r 0 

is dissipative and linear from Xc x to Xc x and since AR : D{Mc) C 
Xc X =1 Xc X is m-dissipative, so is A + Me by Lemma 2.5. Hence there 
is a unique solution (ic, Cx) of equation (6) and we Hnd a unique (x,Xc) € D(A) 
such that (/, fc) + {x, Xc) G A{x, Xc)- □ 


7 Stabilisation via Nonlinear Dynamic Controllers 

The idea of this section is to obtain stability results similar to those for the static 
case, this time in the dynamic controller setup. Our results are based on the idea 
which we employed for the (first) proof of Theorem 5.1 where we took xq G D(A) 
and for x = S{-)xq and some suitable rj G G^([0,1];K) defined 

\ \\x{t)\\l^ + {x{t) I r]PR\{t))^^ . 

Of course, in the dynamic controller scenario we have to add additional terms 
corresponding to the finite dimensional controller Ec as in the preceding section. In 
fact, we assume the following for Ec. 

Assumption 7.1. Assume that A is an m-dissipative operator as in Theorem 6.2 
and further assume that 0 G Mc(0) and there is p > 0 and an orthogonal projection 
n : —>■ on some subspace of such that the following hold. 

1 . 

Re {{{Zc, Uc) I (Xc, '>J-c)))xi,xK^<‘ — ~P (ll^cllxe + Iln'^^clligiVii j , 

{xci'^c} G D(^M(fj^ G ARi^Xc^Uc). 

2. There is a constant c > 0 such that for every [zc, —yc) G AIc{xc, Uc) 

\Vc\^ < c' (llxcllx, + |nuc|^) • (7) 
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3. There are eonstants to,S,c > 0 such that for every mild solution of the non¬ 
linear control system 

l|2^c(io)|lxe ^ (1 ~ l|2^c(0)||x^ + C ||nUc||^2(Q ■ 

Remark 7.2. Note that by time-invariance of the control system the last condition 
then also holds in the form 

lkc(i+ io)|lx, < (l-'5)l|3^c(t)|lx, t>Q. 

Moreover, the fixed time to > 0 may be replaced by kto where k is an arbitrary 
natural number, i.e. in particular we may choose to > 0 as large as we wish. 

Example 7.3 (Collocated case). One particular case which is covered by the pre¬ 
ceding assumption is the following. Let Ec = (Ac, Be, Cc, Dc) be an impedance pas¬ 
sive system with Cc = B* (Hilbert space adjoint with respect to the inner products 
(• I and (• I ■)^Nd on Xc and respectively) collocated, linear and hounded 

and Ac, — Dc he m-dissipative. Further we assume that 

1 . ^c(O) = {0} and there is an equivalent inner product ((• | ■))x on Xc such 
that for some p > 0 and all Xc G D(Ac), Zc G Ac(xc) one has 

Re ((zc I Xc))x, < -P\\xc\\xc ’ (®) 

2. 0 G Dc(0), n : —>■ is an orthogonal projection such that jiiid < 

jllucl (uc G D{Dc), Wc G Dc(uc)) and there is cr > 0 such that for all 
Xc G Xc, Zc G Ac(xc), Uc G and Wc G Dc(uc) one has 

Re (Zc “t" BcUc I Xcfx,. — (CcXc A Wc \ Uc^^Nd er |^c| • (9) 

Then Me = (-Cc -^Dc) satisfies Assumption 7.1. 

Remark 7.4. We give some interpretation for the preceding conditions in the col¬ 
lated input/output case. The Lyapunov condition (8) says that 0 is a globally ex¬ 
ponentially stable equilibrium for the semigroup (Sc(t))t>o associated to Ac. If one 
has a globally exponentially stable minimum at some other point x* G Xc one may 
simply introduce cc"®'" := Xc — x* as new variable to get to the situation as above. 
(Similar, one may choose a nonzero desired equilibrium x(..) Conditions (9) and (7) 
together may he seen as a strict input passivity condition on the controller system 
(after getting rid of the redundant parts of the input which only constitute boundary 
conditions on the system © = (2t, IB, £) In particular if Dc = D). were linear and 
symmetric the second condition would read 

Re (zc + BcUc I Xc)x^ < Re (CcXc + DcUc \ Uc) - a \HdcUc\^ 

for some a > 0 and Hdc Ihe projection on LerDc^. 

We then have the following preliminary, but general result. 

Proposition 7.5. Let © = (21, 25, £) be an impedance passive boundary control and 
observation system and Me ■ D(Mc) Q Xc x ^ Xc x as in Assumption 

7.1. Denote by (S(t))t>o the nonlinear semigroup associated to A as in Theorem 
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6.2. If there is q : X ^ M. such that |( 7 (a;)| < c||x||^ (x G X) and for all solutions 
X € X) n Loo(R+; ^(21)) of X = Six one has q{x) G W^(R+) and 

II^WIIw + ^ : a.e. t > 0, 


then 0 is a globally uniformly exponentially stable equilibrium of {S{t))t>o. 


Proof of Proposition 7.5. Since all the maps S{t) : X x Xc X x Xc are 
continuous, it suffices to consider initial data (x,Xc,o) G D{A). Moreover, we may 
and will assume that to > ^et (xo,Xc,o) G D[A) be arbitrary and 

let (x,Xc)(t) := S(t)(xo, Xc,o) {t > 0) be the unique Lipschitz continuous solution. 
Define 

$(t) := t (||x(t)||^ + ||xc(t)||^J+g(x(t)) + \\xa\\x,ds, t>0 

and note that ^(x,Xc)(t) = (2lx(t), Zc(t)) := M°((x, Xc)(t)) (a.e. t > 0) and <i> is 
Lipschitz continuous and bounded. Then for every t > 2to > 0 we have 


$(t) - $(to) 


= i 


C,Xc)is 


\XxXc 


+ 2sRe (2lx(s) | x(s))j(- +2sRe (zc(s) | Xc(s)) 


Xc 


io 


^ <l(x(s)) + ^ (||Xe(s + to)||^„ - ||Xc(s)||^Jds 


ds 


< 


^0 


lkc(s)||^^ - 2crs |n£x(s)|^ + c ^|*Bx(s)|^ + |n£x(s)|^^ 


- (1 + cc') ||Xc(s)||v 


2 c(l + cc') 


||ne:x(s)||j 


, ds 


I Z /2 ) 

< f (—2as + cc^ + c) |nCx(s)|^ 

Jtn 

c(l+cc) _ 11 (x, X^ ) (s + t Q ) 11 ^ x Xc ) 

) 


< 


< to 


'to 

+ 5a 
c(l + cc') 
6a 

c(l + cc') 


X, Xc)\\L2(^to,2to-,XxXc) ll(^’ ^c)\\L2(t,t+to-,XxXc) 


5a 


||(x,2;c)(to)||jcxXc ■ 


Since behaves as ||(x, Xc)(t)||^ as t —>■ oo, we easily deduce exponential stability 
from this. In fact, for t > to we have 


(x, Xc){t)\\xxXc ~ 


$(t) (7(x(t)) + ||Xc(t)||^_^ 1 2 


- - x 


^ II c||^2(t,i+io} 


^ c\\{x,Xc){t)\\x^x, 


t 


t 


<j^to) + j\\{x,x,){t)f 


XxXc 


< + c') ||(x,x,)(0)||^,^^ + - \\{x,x,){t)\\l^^^ , 
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so that 

||(a;,a;c)(<)||^x^^ < ll(a:, a:c)(0)||5fxXe > max{2to, c'}) (10) 

from where exponential stability with constants M > 1 and w < 0 independent of 
xq follows. From density of D{A) in X x Xc and continuity of S{t) {t > 0) we 
conclude 

||5(t)(a;o,a;c,o)|lxxx, ^ IK^^o,a:c.o)|lxxx„ > (a;o,a:c,o) £ X x t > 0. 

□ 

Lemma 7.6. If & = (21,18, £) is an impedance passive port-Hamiltonian system of 
order N = 1, then there is q € (^^(X;®) with |(;(a;)| < c||a;||.^ {x £ X) such that 
for every solution X € X) D Loo0^+', D{W)) ofx = ‘Xx one has 

II^WIIw + ^'?(a;(t)) < c |('Hx(t))(l)|^ , a.et> 0. 

Proof. Take q{x) = [x \ pP^^^x) ^ as in the proof of Theorem 5.1 and the calcu¬ 
lations made there confirm the assertion. □ 

Thus Proposition 7.5 and Lemma 7.6 together say the following. 

Theorem 7.7. Let G = (21,18, €) be an impedance passive port-Hamiltonian system 
of order N = 1 and Me '■ D{Mc) C Xc x =4 Xc x be as in Assumption 
7.1. Further assume that 

\{'Hx){l)\^ <\^x\^ + \mx\^ , a:£Z)(2l). 

Then the interconnected map A from Theorem 6.2 generates a s.c. eontraction 
semigroup (S{t))t>Q on X X 2Lc with globally exponentially stable equilibrium 0. 


8 Stabilisation of Second Order Systems 

In this section we aim for a generalisation of Theorem 7.7 to the case where 
21a; = Pa (Ha;)" -f Pi(Ha:)' + Po(Hx) 

is a port-Hamiltonian operator of second order (X = 2). Again we assume that 
H £ W)^(0,1)'^^'^ is Lipschitz continuous. For the case of (static and dynamic) 
linear feedback stabilisation, see [2], where exponential stability has been proved 
under the assumption 

|(Hx)(0)|^ -f |(Ha;)(l)|^ -f |(Ha;)'(0)|^ < |25a;|^ -f |ne:a;|^ , x £ P(2l). 

Of course, the proof there used the Gearhart-Greiner-Priiss Theorem, so lacks any 
possible generalisation to the nonlinear scenario. Therefore we aim to apply Propo¬ 
sition 7.5 which amounts to finding a suitable q £ C^(X;R) satisfying the assump¬ 
tions of Proposition 7.5. 

Unfortunately, it is very hard to prove existence of such a functional q without 
any further restrictions on H and the matrices Pq and Pi and, in fact, we did not 
succeed in proving the general result, but had to impose further constraints on the 
matrix-valued function H and the matrices Pq and Pi. 
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Lemma 8.1. Let & = (2t, *8,£) be an impedance passive port-Hamiltonian system 
of order N = 2. Further assume that H',Pq and Pi are small compared to H, i.e. 


2 > +P2'^^’i)(C- 1)1 


Loo(0,l;K‘ix‘i) 


+ ^ + P2"P0 - PlP2^PlP2^)\ 

+ \\\{{P2 "P^rP2 "P^)\\- 


( 11 ) 


Then there is q : X ^ M. such that |( 7 (a;)| < c||x||^ {x € X) and for all solutions 

X € n Loo(R+;-D(2l)) of x = 2lx the function q{x) lies in 

with 

+ ^9(“^(^)) < c(|('Hx)(LO)|^ + \inxy{t,o)f + \inx){t,i)\^) 

for a.e. t > 0. 

Proof. We define the real-valued functional q : X - 


as 


q{x) :=Re (xlrjP^^ J x{^)d^ 


-llH 


3-1 


L 2 


x{m I qPlP2' / , x€X 


( 12 ) 


where the scalar function r] G C''^([0,1];K.) may be chosen suitable at a later point. 
Then for every Lipschitz continuous solution x G fl Loo(R+; L*(2t)) of 

the evolution equation x = 21a; we obtain (omitting the parameter t for brevity and 
employing Lemma 5.4) that 


= Re ( P 2 {nxy' I r]P^^ I x{0d^ ) +Re [x\q j inxf{^)d^ 

L 2 


L 2 


-bRe (^Pi('Hx)'-bPoCHa:) | vPff^ 


+ Re a; I pP^^ / Pi{nxy{0 + PoiPxmd^ 


1/2 


- Re / [Hxyym I vPiP, 


-1 


yCjds. 


L 2 


- Re P. 


3-1 


PiiHxyio + Poynxmd^ I VP1P2 / xiOd^ 


= 2Re {{'Hx)' I qx)j^^ -bRe l^Hx)' \ q' J x{£,)d^^ 

- Re r?(l) ({nxy{l) I [ x{£,)dA - Re ((Px)'(O) | qx)^ 

\ Jo / 

-bRe ( Pi{'Hxy I qP 2 ^ [ x{C)d^] -b Re (a; I 77 / Pf^ PiiHx)'{^)d^ 
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— Re 


x{^)d^ I VP 2 Poi^-x) ] +Re [xlrjP^ Po I {nx){^)d^ 

L 2 


- Re Piinx)' I T]P^^ / + Re Pi(Ha;)'(0) | tjP^ 


-1 


- Re P^'PiCHx) I r]PiR 


+ Re ( P^^Pi{nx){0) I TjPiP- 


L 2 

x{Od^ 


-1 


L 2 

x{C)dI, 


L 2 


L 2 

x{C)d£. 


-Re / P^^Po{nxm\vPiPi 


0 

-1 


L 2 

x{C)dI, 


L 2 


= iivP' - ‘2v'P)x I x)r + [(a;(C) | (?7'H)(C)a:(C))K<i]c 


/ K'^ 


-Re(^x\r]"J x{^)d{^ + r7'(l) Re ^('Ha;)(l) | J x{^)d^j 

-T]{l)Re ({nxYil) \ [ x{^)dA - Re {x \ tjCHxY(0)) 

+ Re (x I r,P,-^P,{nx))^^ - Re (x I A(Ha:)(0))^^ 

-Re x{Y,)d^ \ r]P 2 ^Po{'Hx)^ + Re ^x \ rjP^^ Pq J {'Hx){£,)d£, 


L 2 


+ Re Pi('Ha;)'(0) | t^P^^ / x{£,)d^ 


L 2 


-Re (pi^Pi{nx)\r,PiP^^ [ xiOdA 


+ Re P2'Pi{nx){0) I 77 P 1 K 


-1 


/ L 2 

xiOdY. 


-Re ( P^^PoiUxmivPiP 


-1 


L 2 

x{i)d^ 


L 2 


< {{eU + ry-H' - 2 ri'n + ryRe {P^^PiRYx \ x) 

+ Re (hx I (-77" + P*P-^ 1 ^ + p-^Pov - PiPi^PiPi^v) I xiOd^ 


L 2 


- Re l^P^^Po HxiO I VP 2 Po md^ 
+ c,., ((l + |ry(0)|^)|(?^x)(0)|V|W(0)|^ 
+ Mi)f\inx)ii)f + W{i)f\inx){i)t 


L 2 


(13) 


for every e > 0 and a constant Ce^rj > 0 which may depend on e > 0 and 77, but 
which is independent of x. We now estimate in the following ways. On the one 
hand 

Re fH i I (-77" + P^Pi^v + P^'PoV + P^'PoV - PiP2"PiPOY^v) f x{Od^ 

\ Jo / La 
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< \\nx\\^^ 11 - 77 " + (PSP,-^ + P,-^Po - PlP2'PlP2l4L^iO,l-,K<^^0 

< \\-v" + {P5P2"+Pi^Po - P1P4P1P2" 

and on the other hand 

-Re (p^Po [ {nxmd^ I vP^Po [ xiOd^ 


xiOd^ 




1/2 


L 2 


{Hxmd^ 


L 2 




<\\iP2-'PorP2-^Po4LM:^^xO 
< \\iP2-'PorP2-^Po4LM:^^>^0 ^ ■ 

Therefore, 

tA^'> 

< {[e- 27 j'+ 7 j{n'n-^+ 2 Re{P 4 Pi'H)'H ~4 I 4 

1 - 77 " + iPo*P2" + P 2 " - PiP2"PiP2~")4luotk^>^0 


L 2 


V2 

\\{P2-^PorP2-^Po4L^iO,l;K^x^) 


\\nx\\^j\x\\^^ 


and we only have to find a suitable function 77 such that 

2rj'>e+\\in'n-^+P4Pi)4LUONK^><0 

+ ^ 11 -^" + {PSP2~" + P2~"Po - PlP2"PlP2")4L^iO,l;K^x^) 

+ \UP2-^PorP2-^Po4LUo,i;K^^O- 

In particular for the choice r]{C) = 1 ~ C obtain the condition 

+ ^ ||P*P2-1 + P2-'^"o - PiP4PiP2"\\ 

+ \\\(.P2~"PorP2~"p4 
< 2 . 


The assertion follows. □ 

The interplay of Lemma 8.1 with Proposition 7.5 then implies the following. 

Theorem 8.2. Let & = (21,25, £) be an impedance passive port-Hamiltonian system 
of order N = 2 and Me '■ D{Mc) C Xc x =1 x K'^'^ as in Assumption 7.1. 
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Further assume that 

+ ^ ||(Po*^2-' + Pi^Po - PlP2^PlP2^)\\ 

+ ^||((^2-'^0rP2-'^0)||. (14) 

and 

|CHa;)(0)|^ + \{nx){l)f + |(■Ha;)'(0)|^ < , x G 15(21). 

Then the interconnected map A from Theorem 6.2 generates a s.c. contraction 
semigroup {S{f))t>o on X x with globally exponentially stable equilibrium 0. 


9 Stabilisation of the Euler-Bernoulli Beam 

We investigate how the general result Proposition 7.5 may be used to design ex¬ 
ponentially stabilising controllers for the Euler-Bernoulli Beam equation, i.e. the 
dynamical system 

p{Ou:u{t, C) + {Eiujcc)cdt> C), C e (o, i), t > o. (15) 

The energy of the system is given by 

E{t) = P(C) \ujt{t, 01 " + EI{C) Mt, 01 " dc t > 0. 

To rewrite (15) as a (second-order) port-Hamiltonian system, we set 



and Pq = Pi = 0, so (15) takes the form 

— x(t, 0 = P 2 ^'Hx{t, 0 =: 2ta;(t, 0, C G (0,1), t > 0. 

Note that 

Re {%x I a;)„ = Re [(7^ia;i)'(0*(^2X2)(0 - {HiXi){0*{^ 2 X 2 )'{0]l 

= Re K(0*(1 ^Mc)(C) - cct(0*(l^Mc)c(C)]o 

for all x = (fd() G 75(21). From here, several choices of 25 and £ are possible to 
make © = (21,25, £) an impedance passive port-Hamiltonian boundary control and 
observation system. In that case (and if p~^,EI G LodO, 1) are uniformly positive) 
for any m-monotone ^ =4 the operator A = ^\d{A): P{^) = {x G 75(21) : 

IBx G —(j>{(Lx)} generates a s.c. contraction semigroup on X = L2(0,1)^ (which 
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is a Co-semigroup if ^ G is linear). Lemma 8.1 gives some conditions under 
which the system can be exponentially stabilised, however these conditions are 
rather strong and the proof of Lemma 8.1 does not take into account the additional 
structure of the Euler-Bernoulli beam, in particular those of the matrices Pi. We 
therefore give a result analogous to Lemma 8.1 making use of the Euler-Bernoulli 
beam structure. 

Lemma 9.1. Assume that © = (21, $,£) is an impedance passive second order 
port-Hamiltonian system of the form 



Further assume that Hi G 1 E,^( 0 , 1 ; (i = 1 , 2 ) where d G 2 N is even and 

,\\H[Hf^\\^^} < 1 . 

Then there is q G C^(X;R) with |q(a;)| < c||x||^ (x G X) such that for all solutions 
X G 1 E^(R+; X) n Loo(R+; D{^)) of X = 2 lx one has 

II^WIIw + ^ (l('^^)(0)l^ + I(^i 2 ;i)'( 0 )|^ -b |(H2X2)(1)|^) , a.e. t > 0 . 

Proof. This proof is based on the technique used in [6] for a chain of Euler-Bernoulli 
beams with the particular boundary condition Wt( 0 ) = w^( 0 ) = 0 at the left end. 
Let rj G C^([ 0 , 1 ];]R) be a twice continuously differentiable real function which we 
choose at a later point and define 


q{x) := Re ( xi I 77P W , x = (xi,X2) G X. 


Then for every Lipschitz continuous solution x G lT((o(]R+;X) fl Loo(R+; L*(2l)) of 
the evolution equation x = 2 lx we have (using Lemma 5.4 again) 

S’*"’ 


= Re P *xi,t I ?7 / X2[£.)df, +Re [xi \ r] P X 2 ,ti 0 df, 


= - Re ( {H2X2)” I p / X2{£,)d£, -bRe ( xi I 77 / (Hixi)"{£,)d£, 


= Re {{H2X2)' I 11 x 2 )L, + Re ( {H2X2)' I 77 ' / X2{C)df. 


K‘7/2 


-Re ^77(l)('H2a;2)'(l) | X 2 { 0 d^^ 

-b Re (xi I viHiXi)')^^ - Re (ryxi | ('HiXi)'(O))^^ 

= ix 2 I {{^^■2)' - 277'H^)x2)i, -b ^ [(X2(C) I (77'H2)(C)a:2(C) 
- Re (h2X2 I 77" f a: 2 ( 0 c^'f ) “ Re {H2X2 \ 'n'x 2 )Ln 
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+ Re (?7'(l)('H2a;2)(l) I / - Re (77(1)(H2X2)'(1) | / 2 : 2 ( 0 ^^) 

V 20 /k‘*/2 V Jo /k‘*/2 

- ^ (a^i I iiv'HiY - + i [(a;i(C) | {v'Hi)iOxi{C))K-iMl 

- Re {r]xi I ('Hia;i)'(0))^^ 

< I ii-v' - + £)'H 2 X 2 I X 2 )^^ 

+ ^ {{-v' + + £)'HiXi I - Re (n2X2 \ v" J X2{(,)d^'^ 

+ C, (|,7(l)(772X2r(l) - 7y'(l)(H2X2)(l)|^ + ||77|1^^ \imx,yiO)f 

+ h(0)| |CHiXi)(0)|^ + |? 7 ( 1 )| |('Hia;i)(l)|^) . 

We therefore need to find 77 such that 77(1) = 0 and the following conditions hold 
true for some e > 0. 


v'< 


v'< 


\v\ 


I'X/' "X/”! I 
/to /to 


v"VC 





This gives the assertion of the lemma under the condition on H, using the simple 
choice r]{() = 1 — □ 

Theorem 9.2. Let & = (21,25, £) be an impedance passive port-Hamiltonian system 
of order N = 2 of Euler-Bernoulli type as in Lemma 9.1 and Me '■ D(Mc) Q 
Xc X =4 Xc X as in Assumption 7.1. Further assume that 


|(77x)(0)|" + |(77ixi)'(0)|^ + \{n2X2)il)f < \^xf + inCxl" , X e D{^). 


Then the interconnected map A from Theorem 6.2 generates a s.c. contraction 
semigroup {S(t))t>o on X xX^ with globally uniformly exponential stable equilibrium 


0 . 


Proof. Combine Lemma 9.1 with Proposition 7.5. 


□ 


Acknowledgements 

The author would like to thank Birgit Jacob for arousing his interest in this research 
topic, for fruitful discussions and helpful advice as well as for careful proof-reading. 
He also is grateful to Hans Zwart for making a research visit at the University of 
Twente possible and for several discussions leading to a better understanding of the 
topic and also pointing out a mistake in an earlier version of this article. 


References 

[ 1 ] R.A. Adams: Sobolev Spaces. Pure and Applied Mathematics 65, Academic 
Press, New York 1975. 






PH SYSTEMS WITH NONLINEAR BOUNDARY FEEDBACK 


28 


[2] B. Augner, B. Jacob: Stability and stabilization of infinite-dimensional linear 
port-Hamiltonian systems. Evol. Equ. Control Theory 3 (2014), no. 2, 207-229. 

[3] B. Calvert, K. Gustafson: Multiplicative Perturbation of Nonlinear m-Accretive 
Operators. J. Funct. Anal. 10 (1972), 149-158. 

[4] F. Conrad, J. Leblond, J.-P. Marmorat: Energy decay estimates for a beam 
with nonlinear boundary feedback. Lecture Notes in Control and Inform. Sci. 
147, Springer, Berlin 1990. 

[5] S. Cox, E. Zuazua: The rate at which energy decays in a string damped at one 
end. Indiana Univ. Math. J. 44 (1995), no. 2, 545-573. 

[6] G. Ghen, M.G. Delfour, A.M. Krall, G. Payre: Modeling, stabilization and 
control of serially connected beams. SIAM J. Gontrol Optim. 25 (1987), no. 3, 
526-546. 

[71 K.-J. Engel: Generator property and stability for generalized difference opera¬ 
tors. J. Evol. Equ. 13 (2013), 3II-334. 

[8] K.-J. Engel, R. Nagel: One-Parameter Semigroups for Linear Evolution Equa¬ 
tions. Graduate Texts in Mathematics 194, Springer, New York 2000. 

[9] B. Jacob, H. Zwart: Linear Port-Hamiltonian Systems on Infinite-dimensional 
Spaces. Operator Theory: Advances and Applications 223, Birkhauser, 
Springer Basel 2012. 

[10] T. Kato: Perturbation theory for linear operators. Die Grundlehren der math- 
ematischen Wissenschaften 132, Springer, New York 1966. 

[11] Y. Le Gorrec: Exponential stabilization of boundary controlled port Hamilto¬ 
nian systems with non-linear dynamic control: some preliminary results. (Joint 
work with: H. Zwart, H. Ramirez.) Talk at the workshop Port-Hamiltonian 
Systems: Approximations, Theory and Practice at the Lorentz Center, Leiden 
2014. 

[12] Y. Le Gorrec, H. Zwart, B. Maschke: Dirac structures and boundary control 
systems associated with skew-symmetric differential operators. SIAM J. Control 
Optim. 44 (5) (2005), 1864-1892. 

[13] 1. Miyadera: Nonlinear Semigroups. Translations of Mathematical Monographs 
109, American Mathematical Society, Providence (RI) 1992. 

[14] H. Ramirez, H. Zwart, Y. Le Gorrec: Exponential stability of boundary con¬ 
trolled port Hamiltonian systems with dynamic feedback. IFAG Workshop on 
Control of Sys. Modeled by Part. Diff. Equ. (CPDE), 2013. 

[15] J. Rauch, M. Taylor: Exponential decay of solutions to hyperbolic equations in 
bounded domains. Indiana Univ. Math. J. 24 (1974), 79-86. 

[16] R.E. Showalter: Monotone operators in Banach space and nonlinear partial 
differential equations. Mathematical Surveys and Monographs 49, American 
Mathematical Society, Providence, 1997. 


PH SYSTEMS WITH NONLINEAR BOUNDARY FEEDBACK 


29 


[17] S. Trostorff: A characterization of boundary conditions yielding maximal mono¬ 
tone operators. J. Fund. Anal. 267 (2014), no. 8, 2787-2822. 

[18] J.A. Villegas: A port-Hamiltonian Approach to Distributed Parameter Systems. 
Ph.D. thesis, Universiteit Twente, 2007. 

[19] J.A. Villegas, H. Zwart, Y. Le Gorrec, B. Maschke: Exponential stability of a 
class of boundary control systems. IEEE Trans. Automat. Control 54 (2009), 
no. 1, 142-147. 


